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Abstract. We investigate a fully quantum mechanical spin model for the detection 
of a moving particle. This model, developed in earlier work, is based on a collection of 
spins at fixed locations and in a metastable state, with the particle locally enhancing 
the coupling of the spins to an environment of bosons. The appearance of bosons from 
particular spins signals the presence of the particle at the spin location, and the first 
boson indicates its arrival. The original model used discrete boson modes. Here we 
treat the continuum limit, under the assumption of the Markov property, and calculate 
the arrival-time distribution for a particle to reach a specific region. 
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1. Introduction 

Until recently, in time-of-flight measurements for particles or atoms the quantum nature 
of the center-of-mass motion usually played no role since the particles or atoms were 
very fast. However, the advance of cooling techniques has made it possible to create 
ultracold gases in a trap and produce very slow atoms, e.g., by opening the trap. For 
these low velocities the quantum nature of the center-of-mass motion of an atom can 
have noticeable quantum effects, as the remarkable experiments of Szriftgiser et. al. 
have shown |1 . In the simplest quantum mechanical formulation of a time-of-fiight 
measurement one would create a particle at t = with a localized but extended wave 
function and then ask for the arrival time of the particle at some distant point. Repeating 
this one would get an arrival-time distribution that would depend on the particle's wave 
function. Similarly one might ask for passage or transit times through a region. Such 
questions, and more generally the role of time in quantum mechanics, have attracted 
much interest in recent years [3 El- 
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But how should one measure the arrival time of a particle or atom at some particular 
point and what should the resulting distribution look like? Allcock jl] made an ad hoc 
model of an arrival-time measurement using an imaginary step potential, which leads 
to an 'absorption' of the wave packet; he then identified the absorption rate with the 
arrival-time distribution. In general, this distribution will not be normalized since part 
of the wave packet will be reflected from the imaginary step potential rather than being 
absorbed. Also, part of the wave packet may penetrate the step to some depth before 
being absorbed, thus causing a detection 'delay.' As Allcock noticed, decreasing one 
effect will typically enlarge the other. 

Kijowski proposed physically motivated axioms from which he derived an 'ideal' 
arrival-time distribution for a free quantum particle coming from one direction [5] . The 
resulting distribution agrees with an 'approximate distribution' proposed heuristically 
by Allcock. This distribution has been related [H] to the arrival-time operator of 
Aharonov and Bohm [7j (for more on the latter see |HI and references therein). No 
measurement procedure for the distribution was proposed, and its status, properties and 
generalizations are still being critically discussed in the literature; see e.g. P ITUIITT] . 

Halliwell [T2j employed a detection model based on a single spin coupled to a boson 
bath, a greatly simplified version of a general quantum mechanical detector model that 
was proposed in [T2] and elaborated in [TU Working in one space dimension and 
using Bloch equations he arrived at a Schrodinger equation with an imaginary potential, 
thus giving a basis for AUcock's approach. 

An operational and realistic laser-based approach to the arrival-time problem was 
investigated in [T31[Tni El [11111111201 121] • This approach proposes to measure the arrival- 
time by means of laser induced fluorescence jTHI. The idea is to consider a two-level atom 
with center-of-mass motion, to illuminate some region of space with a laser, and to take 
the detection time of the first fluorescence photon as the arrival time of the atom at the 
(sharp) onset of the laser. In this approach one has to deal with the typical problems of 
delay due to the time needed for pumping and decay of the excited state. There is also 
reflection without detection when the atom is reflected from the laser beam in the ground 
state without emitting a photon. Yet, interesting results could be derived. In the limit 
of a weak laser, there is almost no reflection but a strong delay due to the weak pumping 
to the upper energy level of the two-level system; dealing with this delay by means of 
a deconvolution, one recovers the flux at the position of the onset of the laser from 
the first-photon distribution ^Hj- On the other hand, in the limit of strong pumping, 
reflection becomes dominant and the first-photon distribution is clearly not normalized; 
normalizing it via the operator normalization method of Brunetti and Fredenhagen, 
which preserves the bilinear structure of the distribution j22I; one recovers Kijowski's 
arrival-time distribution from the first-photon distribution |TH]. In this way, Kijowski's 
axiomatic distribution can be related to a particular measuring process. Further, in 
a certain limit it is possible to derive a closed one-channel equation for the ground 
state governing the first-photon distribution ^^l- This equation contains an in general 
complex potential which becomes purely imaginary for zero laser detuning. In this way 
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the fluorescence model makes a connection to AUcock's ad hoc ansatz of an imaginary 
potential. 

In the fluorescence model there is a back-reaction of the measurement on the center- 
of-mass motion of the atom, and this might cause deviations from an ideal distribution. 
It thus seems a good idea to use a measurement procedure that does not interact directly 
with the particle through its internal degrees of freedom, but rather to regard the particle 
only as a catalyst for a transition in a detector or its associated environment. 

Just such a detection model was developed in [31113 El- The model consists of a 
three-dimensional array of D spins (the 'detector') with ferromagnetic interaction. In 
the presence of a homogeneous magnetic field, and for sufficiently low temperature, all 
spins are aligned with the field. Reversing the magnetic field suddenly, such that the 
spins cannot follow the reversal, one can produce a metastable state of this compound 
spin system. The spins are weakly coupled to a bath of bosons. There is a particle 
to be detected and its effect on the collection of spins is to strongly enhance the spin- 
boson coupling when the particle's wave function overlaps that of a detector spin. Thus 
when the particle is close to a spin this spin fiips much faster by virtue of the increased 
coupling to the bath. By means of the ferromagnetic interaction, this in turn triggers 
the subsequent spontaneous fiipping of all spins even in the absence of the particle. 
In this way, the single spin fiip is amplified to a macroscopic event and the associated 
bosons can be measured. The details of the amplification process and the probability of 
false detection due to spontaneous spin fiips were considered in [3 [131 [El ■ The motion 
of the particle, whose presence induces the first spin fiip, was treated classically in the 
calculations. In the following, we will concentrate on the full quantum description of 
this first spin fiip and the quantum mechanical aspects of the particle's motion, and 
comment only briefiy on processes internal to the detector. 

In this paper we investigate this detector model in the limit of continuous boson 
modes, under the condition that the spin-boson interaction satisfies the Markov property 
(see fj43p ) . and use it to determine the arrival-time distribution of a spatially spread- 
out particle. It turns out that one is again led to a Schrodinger equation with an 
imaginary potential and the corresponding arrival-time distribution is similar to that 
of the fiuorescence model. In this detector model there is also a back-reaction on the 
particle of interest. In order to eliminate this back-reaction we discuss the idea of 
decreasing the spin-bath coupling while simultaneously increasing the number of spins. 
It is shown that even in the limit when the spin-bath coupling goes to zero and the 
number of spins to infinity, there remains a back-reaction. 

The plan of the paper is as follows. In Section [2l the detector model is reviewed, 
and in Section [31 the arrival-time distribution obtained from a simplified version of the 
model is calculated by means of standard quantum mechanics. Another approach to 
calculate the arrival-time distribution is presented in Section [H and compared to the 
straight-forward calculation. The advantage of this second approach is that is easily 
extended to the full model, the corresponding calculations shown in [Appendix A[ and 
that it allows to some extent for an analytical treatment of the arrival-time problem. 
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In Section El we discuss the relation of the present detection scheme to the fluorescence 
model. Section IHl deals with the limit of zero coupling and an inflnity of spins, and 
remarks on possible schemes for the optimization of the model and on its application to 
passage-time measurements. 



2. The detector model 

The detector model of jHl El based on the following Hamiltonian. The excited 

state of the j^^ spin is denoted by It)^. and its ground state by ||)^.. Deflne 

^P^ = iT),,(Ti - ii),,ai • (1) 

The Hamiltonian for the detector alone is given by 

H,^, = \ J2 nJ^^a^^ - \ J2 hJj'^a^^ ® ) , (2) 

j j<k 

where huij^^^ is the energy difference between ground state and excited state of the j*^ 

(ik) 

spin, and huuj > is the coupling energy between the spins j and k. 

In addition there is a bath of bosons (e.g. phonons or photons) with free Hamiltonian 

-f^bath = ^ hwidldjt, (3) 

where di is the annihilation operator for a boson with wave vector £. Later a continuum 
limit will be taken. In general, the spins will be coupled to the bath, and there is the 
possibility of spontaneous spin flips due to 

i^spon = Yl ^ (^i'^^'^"' + h-C-) , (4) 

where 

^i^'^ = U),,(TI, ^? = ('^i^'^)' = IT),,ai, (5) 

and the coupling constants and the phases /^"'■' depend on the particular realization 
of the detector and the bath. 

The coupling between the j*^ spin and the bath is assumed to be strongly enhanced 
when the particle is close to this spin. Let the j^^ spin be located in a spatial region 
Qj. The enhancement is taken to be proportional to a sensitivity function x^-^^ (x) which 
vanishes outside Qj, e.g. the characteristic function which is 1 on Qj and zero outside. 
The additional coupling depending on the particle's position is thus 

Hcou, = Yl X^'^ (^) E ^ {9^'^^''''4^-^ + h-c.) , (6) 



with 



9? 

'he 



2 



> 

uU Hami 



2 



tonian is 

H = impart + -f^det + -f^bath + -f^spon + -f^coup , (7) 
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where impart is the free Hamiltonian of the particle, 

impart = pV2m . (8) 

Note that the 'excitation number', i.e., the sum of the number of bosons and the number 
of up-spins, is a conserved quantity. The detection process now starts with the bath in 
its ground state |0) (no bosons present) and all D spins in the excited state |ti • • • Td)- 
As a consequence of the excitation number conservation, it is sufficient to measure the 
state |0) of the bath in order to check whether or not any spin has flipped. For HuJq^ 
only slightly above the energetic threshold set by the ferromagnetic spin-spin coupling, 
and 7^"''' sufficiently small, the probability of a spontaneous spin flip ('false positive') is 
very small |21 CSl EI- But when the particle is close to the j^^ spin, the excited state 
11)^- decays much more quickly, due to the enhanced coupling, 'g^^', of the spin to the 
bath. Then, the ferromagnetic force experienced by its neighbors is strongly reduced, 
and thus these spins can flip rather quickly even in the absence of the particle by means 
of the 7^^; by a kind of 'domino effect', the whole array of spins will eventually flip, 
amplifying the first spin flip to a macroscopic event [3 HSl E] . 



3. The direct approach in the one-spin case 

3.1. A simplified model 

We first consider a simplified model consisting of a particle in one dimension and only 
one spin. This simplification is reasonable if the radius of the region Qj is smaller than 
the distance between spins. (Our assumption of locality of the interaction though is a 
bit stronger than this however, since below in Section for calculational convenience 
we will extend the region Qj to a half-line, i.e., x{x) 0(x).) The vectors x and £ 
are replaced by x and i. Also, we will temporarily neglect ifspon in view of assumption 



7? 



2 



The free 



(i) 
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and accordingly the possibility of spontaneous spin fiips. 
iamiltonian for the particle motion in one dimension is 

HlL = PV2m, (9) 
and the free detector Hamiltonian with only one spin simplifies to 

Hi, = hwo^.. (10) 

The free bath Hamiltonian is given by 

Hltth = ^fi^faYaf^- (11) 

Furthermore, let the spin be located in the interval Id = [0, d] so that 

^cou'p = Xx, (£) E ^ {9i^'^'4^- + h-c.) . (12) 

e 

where the sensitivity function (x) vanishes outside X^. The full Hamiltonian of the 
simplified model is then given by 

trljld rrld , ttI , rrld , rrl,ld /-r q\ 

— -"part "T -"det "r -"bath "t" -"coup- K^'^J 
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This simplified model allows for a direct investigation by means of standard quantum 
mechanics. 

3.2. Energy eigenstates 

To get a first idea of how the present detector model works for an arrival-time 
measurement, we simplify the model in this section a little further by assuming the 
detector to be semi-infinite, extended over the whole positive axis, and take momentarily 

Xi^{x) = e(x), 

where B is Heaviside's step function. Also, we assume for the phases in the coupling 
Hamiltonian fe = throughout this section. The stationary Schrodinger equation with 
energy eigenvalue Ej. for a plane wave coming in from the left, initially no bosons 
present, and the spin in state |t), can be solved piecewise in position space. For x < 0, 
the solution simply reads 



^ki^) = \hr- [e''^ + Mk)e-''^] IT 0) + Mk)e-'''^'^- \l h) , (14) 



where the wave numbers k, ki{k) are fixed by 

h^k"^ fhujQ h^ki{kf hujo , 

= Ek = h hue, 15 

and where |t 0) = |t) |0), || le) = \l) \le)- Note that there is the possibility that 
the particle is reflected from the detector. It may either be reflected after it has 
been detected and a boson of mode i has been created, the coefficient for this event 
being Re{k), or it may even be reflected without being detected, the coefficient being 
Ro{k). The latter will lead to a non-normalized arrival-time distribution. Since this no- 
detection probability is in general momentum dependent the momentum distribution 
of the actually detected part of the wave packet must be expected to differ from that 
of the originally prepared wave packet, hence leading to deviations of the 'measured' 
arrival-time distribution from corresponding 'ideal' quantities. 

For X > 0, the operator H^'^'^ — ]3^/2m is independent of x, because Xi^i^) = ©(a;) 
has been assumed, and it commutes with p^ /2m. The eigenvalues of H^'^'^ —p^ /2m are 
real and denoted by hQ^/2. The corresponding eigenvectors are superpositions of |t 0) 
and II 1^) and denoted by |/x) so that 

{H'''^-py2m)\f^) = ^\^^). (16) 

To obtain an eigenvector of H^'^'^ on x > for the eigenvalue Ek, one has to choose an 
eigenfunction e"^''^*^-'^ of p^ /2m such that 

Ek = {hq^{k)f/2m + hnj2. (17) 

From (|TK|l one has 



Til 

q,{k) = \lk^ + -{u,-%). (18) 



Quantum mechanical detector model for moving, spread-out particles 



7 



Note that qfj,{k) is imaginary if > cuq and 

777 

k'<j (n, - uo) , (19) 

leading to exponential decay. Otherwise qfj.{k) is real. The solution of the stationary 
Schrodinger equation for x > belonging to the eigenvalue can then be written as 

^k{x)-]I^T.'^,{k)e''^^'^''\t^). (20) 

The coefficients q;^(A;), Roik), Re{k) are obtained from the usual matching condition, 
i.e., both 

*^(x) if X < 
*^(x) if x>0 

and its first derivative have to be continuous at x = 0. The eigenvectors can be 
determined numerically. 



- i : ; (21) 



3.3. Detection of a wave packet 

The probabihty of finding the detector spin in state |J,) (and hence the bath in some 
boson state jl^)) at time t is given by integration over the modulus square of the 
respective component of 



Pf-(t)= J2 / dx \{x i (22) 

/oo 
dx \{x T 0|*,)|'^l-Po''^'(t), 
■00 

where the superscript 'disc' distinguishes the discrete model from the continuum limit 
discussed in the next section. As long as no recurrences occur, i.e., no transitions 
I J, 1^) 1-^ It 0), one can regard 

<-(t) = ^^Pf'^t) = -^^o'^^(^)- (23) 

as the probability density for a spin flip (i.e. for a detection) at time t. 
As an example we consider a maximal boson frequency uj^^ and 

ue^u^n/N n=l,---,iV 

ge^ - iG^uJi/N. (24) 

As particle we consider a cesium atom, prepared in the remote past far away from the 
detector such that the corresponding free packet (i.e., in the absence of the detector) 
at t = would be a Gaussian minimal uncertainty packet around x = with Ap and 
average velocity vq. Decomposing this into the eigenstates of H, the wave packet at 
time t is 

/oo 
dkilj{k)^k{x)e-'^''*/^ (25) 
-oo 
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with 

A numerical illustration of wf^^'^{t) for = 40 is given in figure ^ (dots). The numerical 
calculation is time consuming, while in the continuous case with the quantum jump 
approach it is much faster (see next section). 




0.2 
t[[is] 

Figure 1. Dots: spin-flip probability density wf^^'^{t) for an incoming Gaussian wave 
packet of ^ and ^ with Ap = 20/Ltm-i • h and Vq = 1.79m/s; ujq = 2.39 • 10^s-\ 
LU^ = 4.6 • wo, G = 2.782 • lOh-^'^, N = 40. Solid line: wi{t) from ^ for the 
corresponding continuum limit. Up to the time of recurrences, |J, 1^) i— > |t 0) (due to 
the discrete nature of the bath), the discrete and continuum probability densities are 
in good agreement. 



4. Continuum limit and quantum jump approach 
4.1. Basic ideas 

In this section the detector model and its application to arrival-times will be investigated 
in a continuum limit by means of the quantum jump approach [221 • This approach uses 
continuous bath modes as a limit, so that there are no recurrences as in the discrete case. 
It is easily generalized to multiple spins and it is more accessible to analytic treatment. 
The bath modes are eliminated, but in contrast to Bloch equations one can work with 
a (conditional or effective) Hamiltonian and has reduced dimensions. It is based on 
watching for the first appearance of a boson. To do this one would have to observe the 
bath continuously. Since in standard quantum mechanics with the simple von Neumann 
measurement theory this would lead to difficulties associated with the quantum Zeno 
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effect j^|2nil2n|; the quantum jump approach circumvents this by temporally coarse- 
grained observations and a coarse-grained time scale. In the present situation, it reads 
as follows. Instead of continuous observation, one considers repeated instantaneous 
measurements, separated by a time At. For a Markovian system with correlation time 
Tc, one takes At 3> Tc to avoid the quantum Zeno effect, but At much shorter than 
the lifetime of the excited state |Ti • • • Td) in order to obtain a good time resolution. 
Typical numbers for quantum optical models are At ^ 10~^^s . . . 10~^°s. To find no 
boson until t = nAt, no boson must have been found in the first n measurements. The 
probability for this to happen will now be calculated. The detector interval can now 
be finite or semi-infinite. 

Let the complete system (bath, detector, and particle) at to = be prepared in the 
state 

|^o) = |0)|Ti...Td)|^^o), (27) 

where {i/jq) denotes the spatial wave function of the particle. If no boson is found at the 
first measurement then, by the von Neumann-Liiders reduction rule fI7\ I28j. the state 
(up to normalization) right after the measurement is given by projecting with |0)(0|, 

l^fond) = |0)(0| f/(At,0) |0) iTi . . . Td) I^o) , (28) 

where U{t, t') denotes the time evolution operator of the complete system. The 
probability, PQ{At), for no detection is the norm squared of the vector in (j^Hj) . i.e. 

Po(At) =11 |0)(0| f/(At,0)|0)|Ti...TD)|^o) f. (29) 

The state then evolves with f/(2At,At) until the next measurement, and so on. The 
state after the n*'^ consecutive no-boson measurement, j^'cond); i^, up to normalization, 

l^cotd) = |0)(0|f/(nAt,[n-l]At)|0)--- 

■••(0|f/(At,0)|0)|Ti...TD)|V^o). (30) 

The probability, Po{nAt), of finding the bath in the state |0) in all of the first n 
measurements is given by its norm squared, 

Po(r^At) = (v]/^„^,* |vl/^„^,*>. (31) 

Note that (0| [/{uAt, [u — l]At) |0) is an operator in the particle-detector Hilbert space 
which does not rotate |ti • • • Td); by excitation number conservation mentioned after 
(jHl). Thus one can write 

l^coid) = l^cond) = |0) iTl . . . Td) I^Lnd) , (32) 

where t = nAt, and hence 

Poit) ^ {^l""^, I ^^ond) = (^cond | ^cond) (33) 

is the probability that no transition |0) — > i.e. that no detection occurs until the 
time t. The probability for the first detection to occur at next measurement is just given 
by 

Po{t) - Po{t + At) = wi{t)At. (34) 
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The crucial point now is to calculate the 'conditional time evolution' of iV'cond)' i-^-' 
time evolution 'under the condition that no detection occurs', and for this one has to 
evaluate (0| U{uAt, [u - l]At) |0) |ti . . . ]d)- 



4-2. A simplified model 



For greater clarity, the evaluation of (0| f/(z/At, [z/ — l]At) |0) |ti . . . Td) will first be done 
for the simplified model introduced at the beginning of Sectional while the generalization 
to the full model is referred to [Appendix A[ 



We use the interaction picture w.r.t. H^ 



1,1(1 



H 



l,ld 



H^^ and Ujit,t') 



e'n^o ^U{t,t')e 'n^o' Using ()12|) with still discrete, but possibly infinitely many, 
modes a simple calculation gives in second order perturbation theory 

/ 



Uj{uAt, [z/-l]At)|0)|T) = |T) 



11 



v 



dti J dt2 (^(^2)) |^,|2e^(-«--^)(*i-*^) 



(35) 



[u~l]At [u-l]At ^ 

where x{t) = x -\-pt/m is the time evolution of the operator x in the Heisenberg picture 
of the free particle. The phases in the coupling terms have canceled; even if one would 
assume these phases to be dependent on the particle's position, fi{x), this would be 
the case to very good approximation since At is very small and thus x{ti) ~ £(^2) l^ijj - 
Consequently one obtains 

/ 

1- 



Uj{i^At,[u-l]At)\0)\]) = \]) 



\ 



vAt ti 

[u-l]At [i/-l]Ai 

with the correlation function 



K T 



^ g-i(<^<-^o)T 



(36) 



(37) 



To have irreversible decay we go to the continuum limit as follows. At first the 
bath modes are indexed by 'wave numbers' 

i = 27m/Lbath, n = 1, 2, ■ ■ ■ (38) 

and uje is chosen as 

iL!i = c {ue) i, (39) 

so that 

c{ujy ^„ c{ujy 27r 



Au 



c{uj) — ujc'{uj) 



Ai 



C{uj) - UJC'{uj) Lbath ' 



(40) 
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^bath 



11 



(41) 



where r(a;^) does not depend on Lbath- Then one obtains in the continuum hmit by (j4Uj) 



K r 



1 1 



ciu 



ujc'^uj) 



c{ujy 



2 -i(LU-UJo)T 



(42) 



We assume T{uj) to be of such a form that the Markov property holds, i.e. 

k{t) ^0 if r > Tc (43) 

for some small correlation time t^. This is the case, e.g., for r(u;) = F as in quantum 
optics. In the double integral of (IHUj) then only times with ti — 12 < Tc contribute, and 
if Tc is small enough one can write 



Xi, {x (ti)) (x (ts)) ~ Xx, (x {ti)f . 



The double integral then becomes 



At 



dt'x,m' + {i^-i)At)Y 



dr k(t). 



(44) 



(45) 



With At ^ Tc the second integral can be extended to infinity, by the Markov property. 
Putting 

c{uo) - UJoc'{uo) 



A = 2Re 
4hift = 2Im 



dr k,(t) 



UJQ- 



c{ujoy 



|r(^o)r 



(46) 



dr k{ 



T 



one obtains 



UjiuAt,[u-l]At)\0)\]) = IT) 



uAt 



]l-^(A + i53hift) j dtiXx,(x(ti))' 



|t)exp 



[iy~l]At 
vAt 



\ 



^(v4 + i(5shift) 



(47) 



[i/-l]At 

up to higher orders in At. Note that A is a decay rate of the upper spin level; in 
quantum optics A and 5shift correspond to the Einstein coefficient and to a line shift. 
Going back to the Schrodinger picture one then obtains by (jHUjl 



IV^cond) 



e « 



l^o) 



with the 'conditional Hamiltonian 

p"^ h 



Hcond = ^ + 2(^'^hift - iA)X:,^{,Xy 



(48) 



(49) 



Note that this result is independent of the particular choice of At as long as At satisfies 
the above requirements. As a consequence, on a coarse-grained time scale in which At 
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is small, t can be regarded as continuous and iV'cond) obeys a Schrodinger equation with 
a complex potential, 

kLnd) = + o(^shift - iA)x,(xf . (50) 



dt l'^™'^^/ \2m 2 

In this continuous, coarse-grained, time scale, (jH^ yields for the probability density, 
Wi(t), for the first detection 

Mt) = — ^ (51) 
and from (|33|) and (j48p one easily finds 



Wl{t) = ^ (V^cond ^cond - ^^Ld ^cond) 
= A / dxX^XxY\{x l^eond)! • 

Jo 



(52) 



If Xx^i^) is the characteristic function of the interval [0,d] this is just the decay rate of 
the excited state of the detector multiplied by the probability that the particle is inside 
the detector but not yet detected — a very physical result. 

4-3. An example 

As an example we consider the continuum limit of the discrete model of ()24|) . In this 
case one has, with the maximal frequency, 

c{uj) = Co 

= uj^n/N = co27m/Lb^th, n = l,---,N 
^bath = 2ncoN/uj^ 



iG^uJi/N = -iGv^27rcoA 



L 



bath 



Tiu) = <; -iG'^W^^r if u;<u;, 
^ ' ^ else. ^ ' 

In the continuum limit, or Lbath oOj one obtains in the case lj^ > ujq 

|Gp (1 + iu;„jr) e-'(^M-^o)r _ gi^or 



Mr) - ■ ■ ' 



4hift — 2 IG] 

-1 



r2 



f^ln 




-) 









(54) 



and Tc is of the order of uJq . In the integral for wi{t) in (jK^ one has Xx^i^) = ©(a;). 
The resulting wi{t) is plotted in figure ^ for the same wave function and parameters as 
for wf^^^{t) in that figure. Both distributions are in good agreement up to the occurrence 
of recurrences || 1^) i— |t 0) in the discrete case. 
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The agreement is also seen for other values of Ap. If Ax is the width of the wave 
packet in position space and Vq is its average velocity, then the width of the probability 
density for detection is at least of the order of Ax/vq since it takes some time for 
the wave packet to enter the detection region. (Further broadening of the detection 
density arises from the width of the delay of the first spin flip once the particle is 
inside the detector.) Consequently, wave packets with small Ap and thus large Ax 
yield rather broad detection densities. On the other hand, as soon as a significant part 
of the wave function overlaps the detector, in the discrete case the time scale of the 
recurrences is essentially determined by the properties of the detector and the bath 
and by their coupling. In case of wave packets with small Ap long recurrence times 
are needed to obtain a good resolution of the typically broad detection densities. This 
requires a large number of bath modes in the discrete case. We further note that for 
more complicated incident wave packets as, e.g., the coherent superposition of several 
Gaussian wave packets with different mean velocities, the probability density exhibits a 
more complicated structure due to the self-interference of the wave function. 



4.4- The general case 

A procedure analogous to p5|l - (j52|l can be apphed to the three-dimensional model with 
several spins, as explained in [Appendix A| The bosons are allowed to have a direction 
e which varies over the unit sphere. In a continuum limit iV'cond) obeys a Schrodinger 
equation with a complex potential 
Q 

^^di 1^*°°'^^ " ^'"'"^ Wcond) (55) 

where 

i^cond = 1^ + ^{'^shift(x) - iA(x)}. (56) 



A(x) and 5shift(x) are given in ()A.8|) and ()A.9|) . The probability density for the first 
detection is again similar to (jS2I), 

Wl{t) = ^ {^V'cond ^^cond " ^^Ld ^mnd) 

d'xA{^)\{^\^Pl,^)\\ (57) 

which is an average of the position dependent decay rate of the detector, weighted with 
the probability density for the particle to be at position x and yet undetected. 



5. Relation of the present detector model to the fluorescence model 

In the quantum optical fluorescence model [13 HEl 11311311111201 EI] for arrival times one 
considers a two-level atom with ground state |1) = and excited state |2) = (°), which 
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enters a laser illuminated region. In the one- dimensional case one obtains a conditional 
Hamiltonian of the form 

-P' ° ^^^^ ^ (^R) 

^^-'i " 2m + 2 Q{x) -i7 - 2A j ' ^^^^ 

where fl{x) is the (position dependent) Rabi frequency of the laser, A the detuning 
(possibly also position dependent), and 7 the decay constant of the excited level. Note 
that in contrast to the present model this is a two-channel Hamiltonian. The reason for 
this is that the ground state of the quantized photon field is not related to a specific 
internal state of the atom due to the driving by the (classical) laser. In the limit 

» |a E (59) 

where E denotes the kinetic energy of the incident particle, the corresponding 
conditional Schrodinger equation reduces to a one-channel equation for the ground state 
amplitude with the complex potential 

and the excited state can be neglected in this limit p!7]. Physically, condition means 
that the excited state decays very rapidly compared to the time-scales of the pumping 
and the center-of-mass motion. Thus, the first fluorescence photon is emitted, i.e., the 
particle is detected, when and where the excitation takes place. For A = (laser in 
resonance) \^ is a purely imaginary potential, similar as for the detector model outlined 
above; only the physical interpretation of the height of this imaginary potential differs. 
In other words, the one-channel limit of the fluorescence model coincides with the full 
quantum mechanical model from Section |21 when considering the conditional interaction 
for the particle until the first detection. In this way, the fully quantum mechanical 
detector model of Section|2not only justifies the fluorescence model for quantum arrival 
times, at least in the limit of (j59|l . but one can conversely immediately carry over 
the results of the fluorescence model to the detector model. The investigation of the 
fluorescence model has shown that the essential features like reflection and delay |15j . 
and main results like, e.g., linking Kijowski's arrival-time distribution to a particular 
measuring process |lBj, can be obtained from the full two-channel model as well as 
from its one-channel limit. Hence these results immediately carry over to the present 
detector model. Also, the derivation of a complex potential model for particle detection 
from two different physical models, viz. the fluorescence model and the present detector 
model, indicates the importance of the complex potentials and of Kijowski's arrival-time 
distribution, which in turn can be derived from the complex potentials approach. This 
connection is interesting since it can illuminate the physical background of otherwise 
heuristically introduced complex potentials. Differences, however, arise for example in 
applications to passage times since the reset state after a detection is not the same in 
the two models 132 . 
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6. Discussion and extensions 

The model investigated in tliis paper has three ingredients, viz. a particle in whose 
spatial properties one is interested, a 'detector' based on spins, and a bath of bosons, 
originally in the ground state. There is neither a direct measurement on the particle 
of interest nor on the detector but only on the bath, which is checked for bosons. In 
this way one can hope to keep the disturbance of the particle by the measurement to a 
minimum. However, as in the fluorescence model and also seen in ^2] for a simplified 
spin model, also the present full model yields a description by a complex potential 
and thus shows the typical unwanted features: There is a detection delay, due to the 
finite spin decay or flip rate, and there is also necessarily the possibility of reflection 
of the particle by the detector without the detection of a boson, due to the increased 
bath-detector coupling caused by the particle's wave function inside the detector. This 
reflection without boson detection causes a non-detection of the particle so that the 
probability density wi{t) in ()57|1 for the first detection is not normalized. A similar 
effect arises from the transmission of the particle without boson detection. 

In order to reduce the detection delay one may be tempted to increase the spin-bath 
coupling, which mirrors the particle's wave function inside the detector. As a by-product 
this would also decrease transmission without detection. However, the increase of this 
spatially dependent coupling means an increase of the absorbing potential — i^A(x)/2, 
and this will also increase the reflection without boson detection, so much so that in 
the limit of infinite coupling everything is reflected while nothing is detected. The same 
phenomenon occurs in the fluorescence model ^Sl and is a typical feature of complex 
potentials, as already noted by Allcock jl]. 

One can also try to reduce the influence of the spin-bath system on the particle and 
thus the latter's disturbance by decreasing the spin-bath coupling at a space point and 
simultaneously increasing the number of spins located there. This seems natural because 
it is the flip of a single spin which gives rise to the detection, and with a larger number 
of spins this can compensate for the weaker coupling. To investigate this quantitatively 
we consider N spins, later to be taken to oo, in the same volume V and x^''K^) = Xv{^) 
for all j. The coupling constants are taken in the form 



and similarly for 7^ . Further, the ferromagnetic force experienced by the individual 
spin is assumed not to grow with increasing such as for nearest neighbor interaction. 
Then ()A.8|) becomes 




(61) 



N 




)1 r d^]e |r {uo,e)f xvi^y + iTspon (^o,e)|^ 



J (27r)2 N 



I ^2i\^(^o,e)fxv{^? + \rspon{^o,e)f) (62) 
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which is just the decay rate for a single spin in V, with resonance frequency cDq and the 
couphng as for = 1. A similar result holds for 5shift(x), defined in ()A.9|1 . 

Thus, simply increasing the number of spins A^ and scaling the coupling constants 
with \/l/N leaves A and 5shift invariant and thus does not change the dynamics until 
the first detection (spin flip), and in particular does not help to avoid reflection of 
undetected particles. Any other scaling power of A^, however, would not lead to a 
reasonable detector model in the limit N oo since then either A and 5shift would go 
to zero or to cxD. Similar results also hold for the quantum optical fluorescence model. 
It is interesting to note that, although it is the flip of one single spin which triggers the 
detection, it is the totality of all spins located in V which determines the conditional 
time evolution. 

It has been shown in the context of complex potentials, however, that one can deal 
with the delay/transmission- versus-reflection problem by dropping the restriction to 
rectangular potentials [101 EI]. In fact, it is possible to absorb nearly the complete wave 
packet in a very short spatial interval; given a wave packet with a specific energy range, 
an appropriate imaginary potential can be constructed by means of inverse scattering 
techniques. We stress that there is no such a thing as the optimal imaginary potential 
for all wave packets but the construction of the optimized potential requires a priori 
information about energy range of the wave packet under consideration. 

The present detector model is applicable not only to arrival-time measurements, 
but also to more involved tasks like a measurement of passage times. A detailed analysis 
including numerical examples will appear elsewhere [12] • It turns out that a too weak 
spin-bath coupling yields a broad passage-time distribution due to the slow response 
of the detector to the presence of the particle. A too strong spin-bath coupling, 
on the other hand, also yields a broad passage-time distribution due to the strong 
distortion of the wave packet during the measurement process. This is a quantum effect. 
There is, however, an intermediate range for A{x.) yielding rather narrow passage-time 
distributions. Indeed, a rough estimate in shows that for an optimal choice of 
incident wave packet and decay rate A{x.) the precision of the measurement can be 
expected to behave like E~^^^, where E is the energy of the incident particle. For low 
velocities, this means some improvement as compared to the results of models coupling 
the particle continuously or semi-continuously to a clock, where one has £^~^-behavior 
j^E!]. Thus, it appears that the latter behavior of the precision is not due to a 
fundamental limitation related to a kind of time-energy uncertainty relation. 
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Summary 

We have investigated the continuum hmit of a fully quantum mechanical spin-model for 
the detection of a moving particle when the spin-boson interaction satisfies the Markov 
property. In an example with a single spin and 40 boson modes it was shown numerically 
that the continuum limit gave a good approximation to the discrete model up to times 
of revivals. We have derived analytical expressions for the arrival-time distribution. 
The conditional Schrodinger equation governing the particle's time evolution before 
the detection has the same form as the one-channel limit of the fluorescence model, 
which is based on the use of a laser-illuminated region. The quantum spin detector- 
model provides an easier way to obtain this one-channel equation, since no additional 
assumptions or limits are needed. 



Appendix A. The quantum jump approach for several spins 



The continuum limit and quantum jump approach for the full model in Section |21 is 
quite similar until In second order perturbation theory w.r.t. -ffcoup + -f^spon one 

obtains 

/ uAt 

{0| f/,(^A(, \f - 1]A() |0) ITi . . . Td) = ITi " 



wh-E / d, / 

\ ^''^ [u-l]At [u-l]At 



dto 



(A.l) 



where 



D 



(A.2) 

^fc=l k=j+l / 

are modified resonance frequencies arising from the ferromagnetic spin-spin coupling. 
The phases /p'* have canceled similar to the one-spin case since only products of the 
form d£0"^'*a|o"^^ contribute to the second order, and consequently the contributions 
from different spins do not mix. 

Similar to (1421) one can define correlation functions ^2 and K^} in an 

obvious way. Before the continuum limit the bath modes are indexed by the wave 
vectors 

Ui = 1,2, ■■■ (A.3) 



27r 



^bath 



n2 



In analogy to the coupling constants are taken in the form 



-^bath 



(A.4) 
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(A.5) 



with Lj£ = c{u!e)i, e£ = £/i, and 

|r(^-)(a;,,e,)f »|rgLKe,)f. (A.6) 

Again the Markov property is assumed to hold for the correlation functions in the 
continuum limit. The procedure is then analogous to the single spin case, and one 
obtains (j48p with the conditional Hamiltonian 

^2 h. 



Hcond = 1^ + ^{<^shift(x) - iA(x)} 

where ^(x) is given in analogy to ()47j) by 

A(x) =2Re^ / dr{4)(r)x(^)(xf + 4;(r)} 



(A.7) 



(A.8) 



C LU, 



■"0 



UJn C U, 



C UJ, 



W) 



X 



where the df2e integral is taken over the unit sphere and where the contributions 
from 1^}^, have been neglected, due to ()A.6j) . The terms have the familiar form 
of the Einstein coefficients in quantum optics, where there would also be a sum over 
polarizations. 4hift (x) is given by 

3 

Since the Kt^^ term leads to a constant it just gives an overall phase factor and can 
therefore be omitted. 
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